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ABSTRACT 
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We present the results of a Mg II AA2796, 2803 absorption-line survey us- 
ing QSO spectra from the Sloan Digital Sky Survey Early Data Release. Over 
1,300 doublets with rest equivalent widths greater than 0.3 A and redshifts 
0.366 < z < 2.269 were identified and measured. We find that the A2796 
rest equivalent width (W^ 2796 ) distribution is described very well by an ex- 
ponential function dN/dW^ 2796 = ^ e~P" , with N* = 1.187 ± 0.052 and 
W* = 0.702 ±0.017A. Previously reported power law fits drastically over-predict 
^ | the number of strong lines. Extrapolating our exponential fit under-predicts the 

number of Wq 2796 < 0.3 A systems, indicating a transition in dN/dW^ 2796 near 
Wq 279& ~ 0.3 A. A combination of two exponentials reproduces the observed 
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distribution well, suggesting that Mg II absorbers are the superposition of at 
least two physically distinct populations of absorbing clouds. We also derive a 
new redshift parameterization for the number density of Wq 2796 > 0.3 A lines: 
N* = 1.001 ± 0.132 (1 + 2 )°- 226±0 ' 170 and W* = 0.443 ± 0.032 (1 + z )0.634±o.o97 ^. 
We find that the distribution steepens with decreasing redshift, with W* decreas- 
ing from 0.80 ± 0.04 A at z = 1.6 to 0.59 ± 0.02 A at z = 0.7. The incidence 
of moderately strong (0.4 A < W$ 2796 < 2 A) Mg II A2796 lines does not show 
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evidence for evolution with redshift. However, lines stronger than « 2 A show 
a decrease relative to the no-evolution prediction with decreasing redshift for 
z < 1. The evolution is stronger for increasingly stronger lines. Since Wo in 
saturated absorption lines is an indicator of the velocity spread of the absorbing 
clouds, we interpret this as an evolution in the kinematic properties of galaxies 
from moderate to low redshift. Monte Carlo simulations do not reveal any strong 
systematic effects or biases in our results. While more recent SDSS QSO spectra 
offer the opportunity to increase the Mg II absorber sample by another order of 
magnitude, systematic errors in line identification and measurement will begin 
to dominate in the determination of absorber property statistics. 



1 Based on data obtained in the Sloan Digital Sky Survey. 

2 email: dbn@astro.ufl.edu; turnshek@quasar.phyast.pitt.edu; rao@everest.phyast.pitt.edu 
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Subject headings: galaxies: evolution — galaxies: ISM — quasars: absorption 
lines 



1. Introduction 

Intervening QSO absorption lines provide an opportunity to study the evolution of 
galaxies selected via their gas cross-sections, independent of their stellar luminosities, from 
the earliest era of galaxy formation up to the present epoch. The chance alignment of 
galactic gas with the line of sight to a distant QSO allows the unique opportunity to study 
the gaseous phase of galaxies, without the biases involved in luminosity limited studies. In 
particular, low-ion metal lines and strong H I lines can be used to sample the low-ionization 
and neutral gas bound in galactic systems. Surveys for intervening Mg II absorption systems 
are particularly useful, since the Mg II doublet at rest wavelengths AA2796, 2803 allows 
ground based surveys to track them down to relatively low redshift (z = 0.11). Strong Mg II 
systems are good tracers of large columns of neutral hydrogen gas (Rao & Turnshek 2000) 
and are useful probes of the velocity structure of the neutral gas components of galaxies, 
including damped Lja (DLA) systems. Furthermore, unsaturated low-ionization metal lines 
associated with strong Mg II systems can be used to deduce neutral phase metal abundances 
(Pettini et al. 1999; Nestor et al. 2003; Prochaska et al. 2003). 

Many surveys have studied the statistical properties of Mg II absorption systems. 
Lanzetta, Turnshek, & Wolfe (1987) presented the first significant survey and provided the 
benchmark results for the statistics of Mg II absorbers at relatively high redshift (1.25 < 
z < 2.15). They found marginal evidence for evolution in the number density of absorbers, 
but no significant evidence for evolution of the A2796 rest equivalent width (Wq 2796 ) dis- 
tribution, which they fit equally well with an exponential and a power law. Tytler et al. 
(1987) and Sargent, Steidel, & Boksenberg (1988) provided data on Mg II systems at lower 
redshift (0.2 < z < 1.5) and found that the comoving number density of absorbers 3 (dN/dz) 
increases with redshift, consistent within the errors with no evolution. Caulet (1989) com- 
pared the Lanzetta, Turnshek, & Wolfe results to results at lower redshift, finding more 
moderately strong systems and fewer weak systems at (z) = 1.6 as compared to (z) = 0.5. 
The study by Steidel & Sargent (1992, hereafter SS92), which contained 107 doublets with 
rest equivalent widths Wq 2796 > 0.3 A over the redshift range 0.23 < z < 2.06, has served 
as the standard for the statistics of the distribution of Mg II A2796 rest equivalent widths 



3 We use partial derivatives in place of the traditional dN/dz since, as we discuss in §3, the z and Wq 
dependencies of N are not separable. 
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dN/dW^ 2796 , redshift number density dN/dz, and redshift evolution for the past decade. 
Their conclusions included the following: 

1. The W^ 2796 distribution can be fit equally well with either an exponential or power 
law. 

2. dN/dz for systems with W^ 2796 > 0.3 A increases with redshift in a manner consistent 
with no evolution. 

3. The redshift number density dN/dz for the strongest lines does show evidence for 
evolution, with dN/dz decreasing from the no evolution prediction with decreasing 
redshift. 

4. The mean W^ 2796 increases with redshift. 

5. There is no evidence for correlation between the doublet ratio Wq 2796 /Wq 2803 and 
redshift. 

Churchill et al. (1999, hereafter CRCV99), presented a study of 30 "weak" (W* 2796 < 0.3 A) 
systems, and determined dN/dW^ 2796 and dN/dz for the weak extreme of the Mg II Wq 2796 
population. They favored a power law over an exponential for the parameterization of 
dN/dW^ 2796 . 

The Sloan Digital Sky Survey (SDSS; York et al. 2000) provides an opportunity to im- 
prove the statistics of Mg II and other low-ion absorption systems by increasing the number 
of measured systems by orders of magnitude. The statistics and evolution of the absorber 
number density, which is the product of the space density and the absorption cross-section, 
and the distribution of line strengths, which correlate with the number of kinematic subsys- 
tems and absorber velocity dispersion, are important factors in understanding the physical 
nature of Mg II systems and their evolution. Furthermore, the strong connection to DLA and 
sub-DLA systems (Rao & Turnshek 2000; Rao, Turnshek, & Nestor 2004, in preparation) 
makes the parameterization of the Mg II properties key to the study of large H I column 
density systems. Finally, only large surveys provide statistically significant numbers of rare, 
ultra-strong systems for complementary imaging and high-resolution kinematic studies. 

In this paper, we present the results from a Mg II AA2796, 2803 survey using QSO spec- 
tra from the SDSS Early Data Release (EDR; Stoughton et al. 2002). The results of our 
analysis, using only a small fraction of the final SDSS database, is of interest because with 
the large number of Mg II doublets identified (over 1,300 with Wq 2796 > 0.3 A), systematic 
errors already compete with Poissonian errors in the measured statistics. In §2 we describe 
the SDSS EDR data set, our doublet finding algorithm, and tests for systematic biases. In 
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§3 we present the results of our analyses, describe the absorber statistics and their parame- 
terizations, and consider sources of systematic errors. In §4 we discuss our interpretations. 
We present our conclusions in §5. 

2. Analysis 
2.1. The Data 

The SDSS EDR provides 3814 QSO spectra, approximately 3700 of which are of QSOs 
with sufficiently high redshift (0.37 < zqso < 5.30, (zqso) — 1-53) to allow detection of 
intervening Mg II AA2796, 2803 doublets. We analyzed all available spectra, regardless of 
QSO magnitude, as only a small number were too faint to detect the strongest systems in our 
catalog. This data set combines a large number of sightlines with a wide spectral coverage: 
from 3800 A at a resolution of 1800, to 9200 A at a resolution of 2100. This corresponds to a 
Mg II absorption redshift coverage of 0.37 < z abs < 2.27. Typical signal to noise ratios in the 
SDSS EDR QSO spectra are such that the number of detectable systems drops significantly 
for Wq 2796 less than « 0.6 A, though our catalog includes systems down to W A2796 = 0.3 A. 
The strongest system that we found has Wq 2796 = 5.68 A. At all of these strengths, the Mg II 
A2796 line is typically saturated. Thus, no column density information can be gleaned from 
Wq 2796 . However, large Wq 2796 does track high H I columns and exhibits correlation with 
line of sight velocity dispersion and metallicity (see Nestor et al. 2003). Although the EDR 
QSO catalog selection properties are not homogeneous, and therefore not appropriate for 
statistical analyses, the QSO selection criterion should have little affect on the analysis of 
intervening absorption. Effects related to the QSO selection that could impact our results 
are discussed in §3.4. 

2.2. Continuum Fitting 

In order to search the spectra for Mg II doublets, continua were first fitted to the data 
using the following algorithm. Data below the QSO rest-frame Lja emission were excluded 
to avoid the unreliable process of searching for lines in the Lya forest, as were data longward 
of the Mg II AA2796, 2803 emission feature. As our only interest is in detecting and measuring 
the strength of absorption features, we defined a continuum to be the spectrum as it would 
appear in the absence of absorption features, i.e., the true continuum plus broad emission 
lines. An underlying continuum was fitted with a cubic spline, and broad emission and broad 
and narrow absorption features were fitted with one or more Gaussians and subtracted. This 
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process was iterated several times to improve the fits of both the spline and Gaussians. 
Except for a few problematic areas in a small number of spectra, this technique provided 
highly satisfactory results. Note that by avoiding the Lya forest, we do not experience the 
difficulties that arise in trying to define a continuum in that region. The fluxes and errors 
in flux were then normalized by the fitted continua, and these normalized spectra were used 
for subsequent analyses. 



2.3. Mg II Doublet Finding Algorithm 

We searched the continuum-normalized SDSS EDR QSO spectra for Mg II AA2796, 2803 
doublet candidates. All candidates were interactively checked for false detections, satisfac- 
tory continuum fits, blends with absorption lines from other systems, and other special cases. 
We used an optimal extraction method (see Schneider et al. 1993) to measure each Wq\ 

n i 7 \ w _ Ei P ( X i ~ A o) ( l ~ fj) A x m 



v E^ 2 (A*-AoK 
(1 + *W,= V E ,p 2(A ,_ Ao) ^ W 

where P{\ — Ao) represents the line profile, and fa and (jf. the normalized flux and flux error 
per pixel. The sum is performed over an integer number of pixels that cover at least ±3 
characteristic Gaussian widths. Many of the lines found were unresolved. For these, it was 
appropriate to use the line spread function for the optimal extraction profile, and a Gaussian 
with a width corresponding to the resolution generally provided a very satisfactory fit. A 
large proportion of the lines, however, were at least mildly resolved. For these systems, a 
profile obtained by fitting a Gaussian whose width is only constrained to be greater than or 
equal to the unresolved width was appropriate. For the large majority of the mildly resolved 
lines, single Gaussians were satisfactory fits and our method for measuring Wq gave results 
consistent with, and with higher significance than, other methods such as direct integration 
over the line. In rare cases 4%), a single Gaussian was a poor description of the line 
profile, and more complex profiles, such as a displaced double-Gaussian, were employed. 

Identifying Mg II systems requires the detection of the A2796 line and at least one 
additional line, the most convenient being the A2803 doublet partner. We imposed a 5 a 
significance level requirement for all A2796 lines. In addition, we also imposed a 3 a signif- 
icance level requirement for the corresponding A2803 line, in order to ensure identification 
of the doublet. In order to avoid a selection bias in the IF^ 2796 /!^ 2803 doublet ratio {DR), 
which has a maximum value of DR = 2.0, we defined a detection limit for IF A2796 , which 
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we call Wq" 1 , by taking the larger of 5.0 times the error in Wq 2796 or 6.0 times the error in 
Wq 2803 . The errors were computed using the unresolved profile width at the given redshift 
to avoid selection bias in the Gaussian fit-width. Only systems 3,000 km s _1 blueward of 
zqso and redward of Lya emission were accepted. Systems with separations less than 500 
kms" 1 were considered single systems. This value was chosen to match the width of the 
broadest systems found, and is consistent with the window used in the CRCV99 study. We 
also measured Mg I Wq 2852 and Fe II Wq 2600 for confirmed systems, when possible. 



2.4. Monte Carlo Simulations 

We ran Monte Carlo simulations of the absorber catalog to test the efficiency of our 
detection technique and to identify biases and systematic effects. For each detected system, 
a simulated doublet having the same redshift, and similar Wq 2796 , DR, and fit width was 
put into many randomly selected EDR spectra. Each spectrum containing a simulated 
doublet was then run through the entire non-interactive and interactive pipelines, and z and 
Wq 2796 were measured for detected systems. Thus, we had two simulated catalogs: the input 
catalog containing over 9100 simulated doublets, although only ps 4500 of these appear in 
regions of spectra with sufficient signal to noise ratio to meet the detection threshold for 
the input Wq 2796 , and the simulated output catalog, containing ps 4400 doublets recovered 
from the simulation. Measurement error in Wq 2796 and any subtle systematic effects can 
cause systems with Wq 2796 < Wq" 1 to scatter into the output catalog as well as systems 
with Wq 2796 > Wq 171 to scatter out. However, if the input distribution of Wq 2796 does not 
match the true distribution, the ratio of the Pi 4400 recovered systems to the ps 4500 input 
with Wq 2796 > Wq™ is not necessarily an accurate estimate of the biases. Therefore, a trial 
Wq 2796 distribution must be specified according to which lines are chosen from the simulated 
input catalog. The trial distribution must then be adjusted so that the simulated output 
best represents the data. 

As the form of the input distribution is necessarily motivated by the data, further 
discussion of the simulations and their use for estimating systematic errors are discussed in 
§3. 



3. Results 



We were able to identify and measure 1331 Mg II doublets with W^ 2796 > 0.3 A. A small 
number of additional possible doublets were also discovered (< 1% of the sample), but only 
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allowed limits on Wq 2796 to be determined or had questionable identifications and were not 
used in our analyses. 



3.1. W^ 27m Distribution 

Figure 1 shows the distribution of Wq 2796 for our sample. Only those systems with 
Wq 2796 > 0.3 A which we use in our analyses, are shown. The distribution has a smooth tail 
out to W$ 2796 « 4.6 A, with the largest value W^ 2796 = 5.68 A. 

The total redshift path covered by our sample for each value of Wq 2796 is given by: 

/y /» y Nspec 

g(W X2796 , z)dz = J2 M 2796 , z) dz, (3) 

- min J z min j 

where gi (W£ 2796 ,z) = 1 if W$ m (z) < W A2796 and gi (W£ 2796 , z) = otherwise, and the 
redshift limits are defined to be 3,000 km s _1 above Lya emission and 3,000 km s" 1 below 
Mg II emission, or the limits of the data. BAL regions were masked out. The redshift path 
coverage is shown in Figure 2. Figures 1 and 2 were combined to form an unbiased W^ 2796 
distribution, which is shown in Figure 3. The distribution is fit very well by an exponential, 

dN/W™« = ^e"^, (4) 

with the maximum likelihood value W* = 0.702 ± 0.017 A and corresponding N* = 1.187 ± 
0.052. The reduced x 2 comparing the maximum likelihood fit to the binned data is close 
to unity, independent of the choice of bin size. Also shown are data from CRCV99 with 
W X2796 < 0.3 A and SS92 with 0.3 A < W A2796 < 2.9 A. All three data sets have similar 
average absorption redshifts: the CRCV99 data has (z a b s ) = 0.9, and our SDSS EDR and 
the SS92 data have (z a b s ) = 1.1. The SS92 best-fit exponential closely agrees with our SDSS 
results, although our normalization, N*/W*, is slightly (~ la) lower than the SS92 value 
(W* = 0.66 ± 0.11 A, N* = 1.55 ± 0.20). The reason for this offset, as well as the lack of 
scatter as compared to the size of the error bars in the SS92 data, is not clear. We note 
that both our survey and the CRCV99 study consider systems within 500 km s _1 as a single 
system, while the SS92 study uses a 1,000 km s" 1 window. This cannot be the source of the 
normalization offset, however, as it would have the opposite effect, reducing the number of 
systems found by SS92. 

A power law fit of the form ON/dW^ 2796 = CW S with the SS92 values of C = 0.38 
and 5 = 1.65 is shown as a long-dash line in Figure 3. It is a good fit to our data for 0.5 A 
< w x27m < 2.o A. For large values of W^ 2796 , however, the SS92 power law over-predicts 
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ON/dW^ 2796 by almost an order of magnitude. CRCV99 also fit a power law to their binned 
data combined with the SS92 binned data, but excluding the highest SS92 bin. This is shown 
as a dotted line. The combined data sets suggest a transition in dN/dW^ 2796 occurring near 
« 0.3 A. 

We used our simulated catalogs (§2.4) to test for biases in the W^ 2796 distribution. 

We chose lines randomly from the input catalog according to a distribution of the form 

_ w 

dN/dW* 2796 oc e w in , with an initial guess for W* n , until the number of lines recovered were 
equal to that of the actual catalog, which determines N* n . Lines with input Wq 279& < 0.3 A 
were used, although only lines with recovered Wq 2796 > 0.3 A were retained. We determined 
a maximum-likelihood value for W* ut using the recovered Wq 2796 values. We then corrected 
our guess for W* n to minimize \W* ut — W% ata \. The process was repeated several times with 
different seeds of the random number generator to determine variance. 

Except for the weakest systems in our catalog, which were under-predicted, this simu- 
lation was able to match the actual data well. The under-prediction could be a weakness 
of our simulation, as we had few systems with Wq 2796 < 0.3 A with which to model the 
low Wq 2796 end of the distribution. Alternatively, the actual distribution could diverge from 
the simulated exponential for W^ 2796 < 0.3 A, as suggested by the CRCV99 data. Thus, 
though we simulated the entire range W^ 2796 > 0.3 A, we limited the calculation of W* ut to 
W£ 2796 > 0.5 A. Motivated by the CRCV99 results, we then added a second exponential to 
the input distribution while holding fixed the primary distribution. We adjusted the values 
of W* and N* for this second exponential to minimize a x 2 statistic calculated by comparing 
the binned simulated output to the binned data. Thus, we model the input distribution with 
separate "weak" and "strong" components: 

dN/dW^ 2796 = ^ e ~ Wo/w ^ + e ~ Wo/w ^ (5) 

' w: k w: tr • 1 ; 

The resulting best fit values are N* k = 1.71 ± 0.02 and W* k = 0.072 ± 0.001 A, and N* tr = 
0.932±0.011and W* tr = 0.771 ±0.014 A, where the errors are the square root of the variances 
from different choices of random number generator seed. These results are shown in Figure 
4. The "weak" values were very stable under changes in the seed, producing small variances. 
Since the "weak" component was constrained by only a small region of W / A2796 -space in the 
data, the actual uncertainties in the parameters are much larger than the variance. 

Although most of the systems in our sample are at least partially saturated, the dis- 
tribution of DR values (see section §3.5) indicates that the typical degree of saturation 
increases from a mixture at W A2796 < 1 A to virtually all systems being highly saturated at 
Wq 2796 > 2 A. Thus, we note the possibility that if the column density distribution (which is 
not measurable with SDSS data) is described by a power law over the range of line strengths in 
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our sample, curve-of-growth effects could conceivably cause the corresponding dN/dWy 2796 
to deviate from a power law at W^ 2796 > 2 A consistent with our results. Nonetheless, 
the directly measurable dN/dW^ 2796 is parameterized very well by a single exponential for 
systems with 0.3 A < W A2796 < 5 A, and by two exponentials for W^ 27m > 0.02 A. 



3.1.1. Redshift Evolution of N* and W* 

In order to investigate evolution in W*, we split our sample into three redshift bins, 
0.366 < z < 0.871, 0.871 < z < 1.311, and 1.311 < z < 2.269, chosen such that there are 
an equal number of systems in each bin. The three dN / 8Wq 279& distributions are shown in 
Figure 5. Figure 6 shows the resulting W* for each redshift bin (circles) and the Monte Carlo 
input values (squares). The bars represent the variance from the different random number 
seeds only. The curve in Figure 6 is the power law fit described in §3.3. 



3.2. Distribution of Absorption Redshifts 

Figure 7 shows the absorption redshift distribution for our sample. Absorption redshifts 
span the range 0.367 < z < 2.269. The total number of sightlines with sufficient signal to 
noise ratio to detect lines with Wq 2796 > W™ m for several values of W™ m is shown in Figure 
8 as a function of redshift. The conspicuous features at z > 1.5 are due to poor night sky 
subtractions in many of the spectra. The depression near z — 1.1 is due to the dichroic 
(Stoughton et al. 2002.) 

The incidence and variance of lines in an interval of W^ 2796 over a specified redshift 
range are given by: 

dN ^ 1 2 _^ / 1 \ 2 

where the sum is over systems with Wq in the given interval and AZ (Wq) represents the 
path contained in the specified redshift range. Traditionally, dN/dz has been plotted versus 
redshift for lines stronger than a specified W™ m . In Figure 9, we show dN/dz as a function 
of redshift for W^ 27m > 0.3 A, 0.6 A, 1.0 A, 1.5 A, 2.0 A, 2.5 A, 3.0Aand3.5A. Also shown in 
Figure 9 are the no-evolution curves (NECs) for a cosmology with WMAP results (Spergel 
et al. 2003), (Qm, h) = (0.3, 0.7, 0.7), scaled to minimize the x 2 to the binned data. The 
yymin _ 0.3 A, 0.6 A, 1.0 A, and 1.5 A samples have x 2 values that are consistent with no 
evolution. The W™ in = 2.0 A, 2.5 A, and 3.0 A samples are inconsistent with the NECs at 
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> 3(7, while the W™ m = 3.5 A sample is inconsistent at rs 2er. The dotted-boxes in Figure 
9 show the results of the Monte Carlo simulation described in §2.4. 

The large size of our data set allows us to investigate dN/dz not only for distribu- 
tions cumulative in W A2796 , but also for ranges of Wq 2796 . This is potentially important, as 
evolution in the largest W^ 2796 values is not necessarily negligible in cumulative dN/ dz dis- 
tributions. Thus, we repeated the above analysis for the following ranges: 0.3 A < Wq 2796 < 

0.6 A, 0.6 A < w£ 27m < l.o A, l.o A < w£ 27m < 1.5 A, 1.5 A < w£ 2796 < 2.0 A, 2.0 A 

< W X279& < 2.5 A, 2.5 A < W X279& < 3.0 A, 3.0 A < W X279& < 3.5 A, and W X2796 > 3.5 A. 
The results are shown in Figure 10. The 0.6 A - 1.0 A, 1.0 A - 1.5 A, 1.5 A - 2.0 A, 2.0 A - 
2.5 A and 2.5 A - 3.0 A samples have x 2 values that are consistent with no evolution. The 
NEC for the 0.3 A - 0.6 A sample is ruled out at « 2.5a and for the 3.0 A - 3.5 A and > 3.5 A 
samples at « 2.0a. 

Since the NEC normalization is a free parameter, a plot cumulative in redshift comparing 
dN/dz from the data to the NEC is more instructive. These plots, shown in Figures 11 and 
12, highlight the skew of the dN/dz curves from the NEC prediction that is not necessarily 
manifested in the x 2 analysis. The NECs over-predict dN/dz at low redshift for small 
W™ m and under-predict dN/dz at low redshift for large W™ m . The transition occurs for 
the cumulative W^ 2796 plots between W™ m = 0.6 A and 1.0 A and the effect is stronger for 
increasingly larger values of W™ m . The plots using ranges of Wq 2796 show that the transition 
from over- to under- predicting dN/ dz at lower redshift occurs around 1 A, and significant 
detection of evolution is seen in systems with Wq 2796 > 2.0 A. The evolution signal is strong 
4a) for lines with Wq 2796 > 3.5 A. Although these plots do give a clearer indication of the 
deviation from the NECs, it is difficult to ascribe a K-S probability to the curves in Figures 
11 and 12 because systems contribute to dN/dz in a non-uniform manner, i.e., inversely 
proportional to g(Wo, z). 

The dN/dz point in the lowest redshift bin for the Wq 2796 > 0.3 A sample in Figure 9, 
and for the 0.3 A < W X2796 < 0.6 A sample in Figure 10, lie well above the NEC. The increase 
is greater in magnitude but smaller in significance (~ 2a versus ~ 3a) for the non-cumulative 
sample. The Monte Carlo results lessen the significance in the cumulative sample, but not in 
the non-cumulative sample. Thus, the weakest lines in our study show evolution in the sense 
that their incidence increases with decreasing redshift. However, this may be an artifact of 
the sharp cutoff in our sample at 0.3 A. A sample with Wq 2796 > 0.4 A is consistent with the 
NEC, and dN/dz for 0.4 A < W^ 279& < 0.6 A is within la of the NEC at all redshifts. We 
cannot exclude the possibility, however, that the increase is real. 

A comparison of the dN/dz values from our survey and SS92 for the full redshift ranges 
are shown in Table I. Our Wq 2796 > 0.6 A and 1.0 A results are consistent with SS92. 
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However, we find a smaller value of dN/dz for W A2796 > 0.3 A (although this difference is 
confined to z > 1.3.) 

3.3. Joint H/ A2796 -Redshift Distribution 

Conventionally, the W^-redshift distribution of absorption lines, d 2 N/dz dW , has been 
parameterized by a combination of a power law in redshift and exponential in W (e.g., 
Murdoch et al. 1986; Lanzetta, Turnshek & Wolfe 1987; Weymann et al. 1998): 

N* w Q 

dN/dWo = 7T^e-w*, dN/dz = N (1 + z)\ (7) 
W* 

This parameterization is convenient for the determination of best-fit values for the parameters 
because of the separate W$ and z dependences. However, our data reveal that a single 
power law is not a good fit to dN/dz for any given range of Wq 279& . Also, we have shown 
that W* depends on redshift for Mg II A2796 lines (i.e., dN/dz for different ranges of 
Wq 2796 varies differently with redshift.) We do not wish to abandon the exponential form, 
however, as we have also shown that it is a very good parameterization of the data at all 
redshifts for the range of Wq 2796 covered by our sample. Thus, we retain the general form 
dN/dWy 2796 = , but allow both N* and W* to vary with redshift as a power law 

in (1 + z). The maximum likelihood results are N* = 1.001 ± 0.132 (1 + z )°- 2 26±o.i7o and 
W* = 0.443 ± 0.032(1 + z )°-634±o.097 ^ The va i ues are highly correlated, and the errors 
include the effects of the correlations. See the Appendix for details of our parameterization 
and a discussion of the errors. 

We repeated the analyses of §3.2, comparing the dN/dz curves to those determined 
using our redshift parameterization of N* and W*. The results, shown in Figures 13 and 
14, indicate that our parameterization is indeed a good description of all of the data. The 
evolution is well described by a steepening of the Wq 2796 distribution with a normalization 
that changes so as to keep dN/dz constant for the smaller Wq 2796 values. Extrapolating 
our parameterization to lower redshift, we predict dN/dz = 0.28 ± 0.03 for Wq 2796 > 0.6 A 
and (z) = 0.06 and dN/dz = 0.11 ± 0.02 for W 7 ^ 2796 > 1.0 A and (z) = 0.04, in comparison 
to Churchill's (2001) values of dN/dz = 0.22+°;J 2 and dN/dz = 0.161^, respectively. 
However, our parameterization is likely inappropriate for describing H 7 ^ 2796 < 0.3 A lines, 
since not including the apparently two-component nature of the Wq 2796 distribution leads to 
underestimates of the number density of weak lines (Figure 3). 

Table 2 summarizes the various parameterizations of the data. 
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3.4. Systematic Errors 

Our Monte Carlo simulations did not uncover any significant systematic errors. All of 
the Monte Carlo dN/dz values are consistent with those from the data. However, since we 
used the data to model the simulated lines, if there are lines with characteristics that make 
them underrepresented in the data, they would also be underrepresented in the simulated 
catalogs. We ran simulations to address this issue as we developed our line-finding algorithm. 
Although it is possible that we miss certain types of lines, such as weak kinematic outliers 
or lines with abnormally broad or exotic profiles, for example, these effects are likely to be 
small. 

A potentially more serious source of systematic error may arise in the regions of poorly 
subtracted night sky lines seen in some of the spectra. The error arrays are not always accu- 
rate in these regions, and doublets falling between sky lines can be confused with residuals 
from the poor subtraction. As the errors are larger (though not necessarily accurate) in the 
night sky regions, the values of Wq" 1 are large as well. Thus, only the strongest lines would 
be affected. In principle, our simulations should account for these effects. However, since 
non-Gaussian profiles are preferentially found among the stronger lines, their simulation is 
somewhat less reliable. Also, their numbers are much smaller, providing fewer lines to serve 
as models in the simulation, and less overall significance. The largest Wq 2796 ranges were, 
in fact, too sparse for meaningful simulations. Errors of this type would be manifest in the 
largest Wq 2796 ranges and for redshifts z ~ 1.6 — 2.0. 

Finally, in QSO absorption-line studies it is of interest to asses if biases due to the pres- 
ence of the absorber affecting the magnitude and color of the background QSO are present. 
The most often discussed effect is that of a dimming and reddening of the background 
QSO due to dust in the absorber. However, the presence of an absorbing galaxy could also 
contribute light and/or have a lensing effect, causing the QSO to appear brighter. These 
competing effects are investigated, using the absorbers and simulations from this work, in 
Menard, Nestor, & Turnshek (2004). They find that QSOs with Mg II absorbers tend to be 
brighter in red passbands and fainter in blue passbands, indicating a combination of bright- 
ening and reddening. The effects are stronger for increasingly stronger Wq 2796 systems. The 
average effects are mostly within ±0.2 magnitudes. As the SDSS EDR QSO catalog selection 
properties were not necessarily homogeneous, it is difficult to quantify the effect these results 
may have on our data. If these effects are non-negligible, they could affect the measured 
WA2796 distribution and dN/dz. The biases are likely to be small, however, since it would be 
peculiar if they conspired to drive dN/dW^ 2796 to more closely resemble an exponential, and 
there is no detectable deviation from an exponential, consistent with the expected biases, 
seen in Figure 3. Since dust increases with decreasing redshift (consistent with the findings 
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of Nestor et al. 2003), it is worth considering if the evolution in dN/dz seen for the strongest 
systems is an affect of increasing dust. However, Menard et al. find no such increase in bias 
with decreasing redshift. Also, Ellison et al. (2004) compare Mg II dN/dz determined from 
radio-selected CORALS survey spectra to values obtained with optical surveys (including 
this work), and find good agreement (though they do not investigate the largest Wq 2796 
ranges). Therefore, the steepening of the Wq 2796 distribution due to a disappearance of the 
largest Wq 2796 lines at low redshift and the evolution in dN/dz for strong lines appears to 
be real. 

While more recent SDSS QSO spectra offer the opportunity to increase the Mg II 
absorber sample by another order of magnitude, systematic errors in line identification and 
measurement will begin to dominate in the determination of absorber property statistics. 



3.5. Mg II Doublet Ratio 

Mg II doublet ratios span the range from DR = 2.0 for completely unsaturated systems 
to DR = 1.0 for completely saturated systems. Figure 15 shows the DR distribution for our 
sample. The doublets are, for the most part, saturated. Ninety percent of doublets with a 
measured DR are within 3<r of DR = 2.0, and 92% have DR — 1.0 > a DR . Figure 16 shows 
DR as a function of redshift. There is no detectable evolution in the DR distribution. In 
fact, the DR distribution is remarkably consistent over the three redshift ranges of Figures 
5 and 6. 



3.6. Mg II Velocity Dispersions 

For absorption lines that are at least partially saturated, such as those that dominate 
our sample, Wq 2796 is primarily a measure of the number of kinematic subcomponents along 
the line of sight (Petitjean & Bergeron 1990) and, to a lesser extent, projected velocity 
dispersion, <J ve i tP (see, for example, Churchill et al. 2000, hereafter C2000). We can directly 
extract information about u ve i, P for individual systems by deconvolving the fitted profile from 
the line spread function. However, as most of the A2796 lines are only mildly resolved, this 
approach is not sensitive for all but the most resolved systems. Nonetheless, we investigated 
the evolution of the a ve i ;p distribution over the three redshift ranges of Figures 5 and 6. In 
all three redshift ranges, ~ 80% of the systems have measured values of a ve i, P < 90 km s -1 . 
Comparing just systems with o ve \ jP > 90 km s _1 , we find that a ve i ;P tends to be smaller in 
the lowest redshift bin (see Figure 17). However, the middle- and lower-redshift bin velocity 
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distributions only differ at the la level. 

3.7. Fe II A2600 and Mg I A2852 

Figure 18 shows the distribution of Fe II Wq 2600 for our sample. The error-weighted 
mean is (Wq 2796 /Wq 2600 ) = 1.42 and the distribution has a smooth tail extending out to 
large values, though ratios above ~ 4 are dominated by values with significance less than 3c. 
Figure 19 shows the distribution of Mg I 1F A2852 for our sample. The error-weighted mean 
is (lF A2796 /lFo 2852 ) = 4.14 and the distribution has a smooth tail extending out to large 
values, though ratios above ~ 8 are dominated by values with significance less than 3a. 

4. Discussion 

The determination of d 2 N/dz cW A2796 provides insight into the nature of the systems 
comprising the population of Mg II absorbers. Specifically, the inability of a single functional 
form to describe dN/ dW^ 2796 suggests that multiple physical populations contribute to the 
absorption. Additionally, dN/ <9W A2796 can be used to infer average absorption cross sections 
for a given Wq 2796 regime. The W / ( ^ 2796 -dependent evolution of dN/dz holds further clues 
related to the nature of Mg II absorbers. These issues are discussed below. We also note 
that a significant fraction (« 36%) of Mg II systems with W / " A2796 > 0.5 A and 1F A2600 > 0.5 A 
are DLAs. The implications of the evolution in the Mg II dN/dz to Vt DLA will be discussed 
in Rao, Turnshek, & Nestor (2004). 

4.1. The Nature of the Absorbers 

The distribution of W / " A2796 for absorbers with W A2796 > 0.3 A is fit very well by a single 
exponential. However, extrapolating this exponential to smaller values of Wq 2796 under- 
predicts the incidence of lines (Figure 3), motivating a two-component description of the 
form of Equation 5. This is the first clear indication of such a transition and raises the 
question of whether the ensemble of clouds comprising Mg II absorbers are the result of two 
physically distinct populations. 

In a series of papers using HIRES/Keck data (Charlton & Churchill 1998; C2000; 
Churchill & Vogt 2001; Churchill, Vogt, & Charlton 2003), Churchill and collaborators have 
investigated the kinematic structure of Mg II absorbers at 0.4 < z < 1.2. For intermedi- 
ate/strong absorbers {Wq 2796 > 0.3 A), they show that the absorption systems are composed 
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of multiple kinematic subsystems, usually containing a dominant subsystem (ensemble of 
clouds) with a velocity width m 10 — 15kms -1 and a corresponding H I column density 
which is optically thick at the Lyman limit (i.e., Nhi >3x 10 17 atoms cm -2 ). Additionally, 
a number of weaker, often unresolved (at ~6 km s -1 resolution) subsystems (clouds) with a 
large spread of velocity separations (up to ~ 400 km s -1 ) from the systemic velocity of the 
absorber are usually present (see Churchill & Vogt 2001, Figure 7.) These weaker subsys- 
tems typically contribute only 10%-20% of the rest equivalent width. The strongest Mg II 
absorbers (Wq 2796 > 2 A) often contain more than one of the dominant subsystems and have 
equivalent widths dominated by saturated features. 

CRCV99 presented a HIRES/Keck study of weak (1V A2796 < 0.3 A) Mg II absorption 
systems. The properties of the kinematic subsystems comprising weak Mg II absorbers are 
similar to the weak subsystems of intermediate/strong Mg II absorbers. The individual 
clouds have W^ 2796 < 0.15 A and sub-Lyman limit H I columns. 

Although most of the Mg II A2796 lines in our survey are at least partially saturated, 
the degree of saturation is correlated with W A2796 . For example, the DR values for systems 
with 0.3 A < W£ 2796 < 1 A indicate a mix in degree of saturation. For lines in this regime, 
dN/dWy 2796 is sensitive to the distribution of Mg II column densities, dN/dN Mg ii- In fact, 
Churchill & Vogt (2001) and Churchill, Vogt, & Charlton (2003) find almost identical power 
law slopes for ON/dW^ 2796 and dN/dN Mg n in their 0.3 A < W^ 2796 < 1.5 A sample. Systems 
with Wq 2796 > 2 A, however, are almost always highly saturated. Also, a general distinction 
between weak (Wq 2796 < 0.3 A) and intermediate/strong (W{f 2796 > 0.3 A) Mg II absorbers 
is the absence or presence of a "dominant" subsystem, though the differences are neither 
discrete nor without exception. Additionally, Rigby, Charlton, & Churchill (2002) note a 
substantial excess of Mg II absorbers that are comprised of a single kinematic component, 
which are almost exclusively weak absorbers. Therefore, given the form of dN/dW^ 2796 , it 
is appropriate to consider a picture in which the sum of the weaker kinematic subsystems 
(clouds) in an absorber are described by N^ k and W* k , and the stronger (dominant ensem- 
bles) subsystems by N* tr and W* tr . If this multi-population explanation is indeed correct, the 
nature of these populations needs to be explained. First, we consider the weak subsystems. 

It has been demonstrated that absorption line systems can be associated with galaxies 
out to large galactocentric distances (~ 100/i -1 kpc for C IV and ~ 200/i -1 kpc for Lya 
absorbers, for example; Chen, Lanzetta, & Webb 2001; Chen et al. 2001). Though there 
is evidence for Lya absorbers in voids (Penton, Stocke, & Shull 2002), these appear to be 
limited to the weakest systems. Also, Rigby, Charlton, & Churchill (2002) conclude that a 
significant fraction of relatively strong (10 15 ' 8 cm -2 < N H i < 10 16 ' 8 cm -2 ) Lya forest lines 
are associated with single-cloud weak Mg II absorbers, which are also C IV absorbers. 
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Therefore, we consider whether single-cloud/Lya forest absorbers at relatively large 
galactocentric distances can contribute significantly to our "weak" component of dN/dWQ 2796 . 
CRCV99 claim to be 91% complete down to W^ 2796 = 0.03 A. Thus, using; the measurements 
of single-cloud systems from CRCV99 with W^ 2796 > 0.03 A, we find W* c = 0.063 ± 0.015 A. 
Rigby, Charlton, & Churchill (2002) give a dN/ dz value for the single-cloud systems, which 
can be used to normalize dN sc / 0Wq 279& . While our comparison is approximate in that the 
redshift range covered by the C2000 sample, 0.4 < z < 1.4, is different than that of our 
sample, the result shown in Figure 20 demonstrates that the single-cloud/Lya forest sys- 
tems should contribute significantly to the upturn in dN/dW^ 2796 , especially for the weakest 
systems. 

Although the Mg II absorber galaxy population has been shown to span a range of 
galaxy colors and types, Steidel, Dickinson, & Persson (1994) describe the "average" galaxy 
associated with an intermediate/strong Mg II absorber galaxy as one which is consistent 
with a typical 0.7L* Sb galaxy. Furthermore, Steidel at al. (2002, also see Ellison, Mallen- 
Ornelas, & Sawicki 2003) compare the galaxy rotation curves to the absorption kinematics 
for a sample of high inclination spiral Mg II absorption galaxies. They find extended-disk 
rotation dominant for the absorption kinematics, though a simple disk is unable to explain the 
range of velocities, consistent with a disk plus halo-cloud picture. However, there are several 
counter evidences to a rotating disk description for the dominant subsystems. The most 
important is that imaging often fails to reveal a disk galaxy in proximity to the absorption 
line of sight. DLA galaxies, which are a subset of strong Mg II absorbing galaxies, are 
not dominated by classic spirals. Of the 14 identified DLA galaxies summarized in Rao 
et al. (2003), only six are spirals. Furthermore, systems that have been studied without 
the discovery of the absorbing galaxy despite deep imaging rule out at least bright spirals 
(although disky-LSBs may still contribute.) Bright galaxies near QSO sightlines usually show 
Mg II absorption, but not all strong Mg II absorbers have a nearby spiral or bright galaxy. 
Thus, though it is clear that rotating disks do contribute to the population comprising the 
dominant kinematic subsystem(s) of Wq 2796 > 0.3 A absorbers, they can only account for a 
(perhaps small) fraction of the total population. 

The detailed nature of the remaining contribution is yet unclear, but this dichotomy 
supports the idea of multiple populations. Bond et al. (2001), for example, show that the 
kinematics of many strong (defined as Wq 2796 > 1.8 A in their sample) Mg II absorbers have 
kinematic structure that is highly suggestive of superwinds/superbubbles. They also show 
that DLAs exhibit low-ion kinematic profiles distinct from the superwind-like absorbers, 
which tend to be sub-DLA. Furthermore, imaging studies (e.g., Le Brun et al. 1997; Rao 
et al. 2003) have found LSB, dwarf, and groups/interacting systems, in addition to some 
spirals, responsible for strong Mg II absorption systems (though the samples are typically 
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selected to be DLA sightlines). Sightlines through groups and interacting systems are likely 
to sample large velocity spreads. Apparently, their contribution to the total cross section of 
strong Mg II systems is non-negligible. 

Thus, though a two-component description is sufficient to describe <9iV/cW A2796 within 
the current data, it is likely that more than two physically distinct populations contribute 
significantly to the phases comprising the absorbing clouds. By observing a large number of 
galaxy types that sample a range of W^ 2796 and impact parameter, it would be possible to 
ascertain what populations contribute most directly to the absorber population. 



4.2. The Absorber Cross-Sections 

The relative mean absorber cross-section, a(z), for ranges of Wq 2796 can be determined 
directly from ON/OWq 2796 . Additionally, knowledge of dN/dz and the galaxy luminosity 
function (LF) allows the approximation of the normalization of the mean absorber cross- 
section, under the assumption that the absorbers arise in galaxies, since dN/dz is the prod- 
uct of a(z) and the number density of absorbers, n(z). However, the process is inherently 
uncertain since: (a) the redshift dependence of the galaxy LF is not well determined; (b) 
although a scaling relation a(L) oc (L/L*) 13 has traditionally been assumed, the determina- 
tion of P is uncertain and likely not universal for all morphological types or redshifts, if it 
is even appropriate; and (c) depending on the value of (5 and the LF faint-end slope, the 
minimum absorber galaxy luminosity (which also is likely to be dependent on morphology 
and redshift) can have large effects in the statistically-determined cross sections. 

Nonetheless, an approximate normalization is of physical interest. Thus, we define R e ff 
to be the effective projected radius for absorption such that <r(L) = nRl^(L). Note that 
the galactocentric radius within which the absorption may occur (i.e., the range of impact 
parameter) may be very different than R e /f- Steidel, Dickinson, & Persson (1994) find 
absorbing galaxies as faint as L ~ 0.05 L*, and a scaling law for the gaseous cross-section 
R e ff(L)/R*ff = (L/L*) 2 for both optical (B-band) and near-IR (K-band) luminosity. They 
claim the relation is much tighter for the near-IR. Thus, we use the redshift parameterization 
of the K-band LF from the MUNICS data set (Drory et al. 2003) and L min = 0.05 L*, though 
we also test L min = 0.25 L* and L min = 0.001 L* to investigate the effects of using high/low 
values of L rnin . We adopt the Dickinson & Steidel (1996) revised scaling law for Mg II 
systems with W^ 2796 > 0.3 A, namely R eff (L K )/R* eff = (L K /L* K ) - 15 . 

The results for the different W™ m values from §3.2 are shown in Figure 21. Despite the 
limitations of this approach, approximate values for R* e *f are apparent: ~ 60 — 100 kpc for 
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W^ 2796 > 0.3 A, « 30 - 60 kpc for W^ 27m > 1.0 A, « 15 - 30 /i^ 1 kpc for W 7 ^ 2796 > 2.0 A 
and < 10 /i^ 1 kpc for W 7 ^ 2796 > 3.5 A. Guillemin & Bergeron (1997) find = 65 /i^ 1 kpc 
using the upper envelope of the distribution of observed impact parameters for a sample with 
0.3 A < W^ 2796 < 4.7 A and (W$ 279e ) = 1.3 A. CRCV99 find R* eff = 110 h m l kpc for their 
w \2796 > 0.02 A sample, but use a LF normalization that is ~ 3.5 times our adopted values. 
Chen, Lanzetta, & Webb (2001) find R* e ^ = 137 kpc for C IV absorbers, which may 
sample larger galactocentric radii. The R* e ^ values are averages; it is expected that there are 
large galaxy to galaxy variations. Direct comparison of R*^ to the distribution of impact 
parameters would require knowledge of the gas geometry and the (likely W A2796 -dependent) 
covering factor. For example, small regions that allow large Wq 2796 absorption may be found 
at galactocentric radii r 3> R e ff if the covering factor for such absorption is / <SC 1. 

We also note that for the single-cloud Mg II systems (see §4.1), Rigby, Charlton, & 
Churchill (2002) derive dN/dz = 1.1 ± 0.06. They claim that these systems are 25%-100% 
of the Lya forest with 10 15 8 cm -2 < N(R I) < 10 16 ' 8 cm~ 2 . Under the assumption that they 
can be associated with galaxies, their incidence corresponds to R*^ ~ 80 — 120 kpc in 
the above approximation. This can be compared to the claim of Chen et al. (2001) that all 
iV(H I) > 10 14 cm~ 2 Lya forest lines arise within the characteristic radius R* ~ 250 kpc 
of luminous galaxies. 

4.3. Nature of the Evolution 

4.3.1. W^ 2796 < 0.3 A Systems 

The Lya forest dN/dz is known to decrease strongly over the redshift interval 4 > 
z > 1.6, but little evolution is present at lower redshifts (see Weymann 1998, and references 
therein). In §4.1, we speculated that a portion of the enriched Lya forest contributes to the 
weak subsystems found in Mg II absorbers. If this is indeed the case, then this evolution 
may manifest itself in dN/dz for the weakest (W 7 ^ 2796 < 0.3 A) Mg II systems at high z. 

4.3.2. The Lack of Measured Evolution in Intermediate/ Strong (0.3 A < W^ 2796 <2A) 

Systems 

Our sample covers the redshift interval z = 2.267 to z — 0.366 which corresponds to 
6.6 Gyrs, or about half the age of the universe. As long assumed and recently demonstrated 
with the Hubble deep fields (see Ferguson, Dickinson, & Williams 2000 for a review), galaxies 
were quite different at z — 2.3 than at z — 0.4. It is therefore noteworthy that, even with 
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the large absorber sample presented here, relatively little evolution is detected in the bulk 
of Mg II absorption systems. For example, it is believed that the global star formation rate 
peaked near z = 1 and has declined by a factor of ~ 2 by z = 0.5 and by a factor of ~ 10 
by z = (Hopkins 2004 and references therein.) Contrastingly, the total cross-section for 
absorption n(z)a(z) of our W™ m = 0.6 A and 1.0 A samples evolves < 10% from z = 1.9 
to z — 0.5 (corresponding to an interval of 5 Gyrs). The corresponding R e ff values do, 
however, indicate that much of this cross section is at large galactocentric radius, extending 
well beyond stellar galactic radii. 

Two immediate conclusions can be drawn from the lack of evolution in the total cross 
sections. First, a large majority of the structures responsible for the bulk of the absorp- 
tion cross-section were in place by z ~ 2. Second, either the time scales governing the 
dynamics of the structures is greater than several Gyrs, or some process(es) regulate the 
production/destruction of the structures such that a nearly steady state is reached. 

Mo & Miralda-Escude (1996) explore two-phase models for gaseous galactic halos in 
which "cold phase" clouds condense from the extended "hot phase" ionized halo at a cooling 
radius r c , and undergo infall at v ~ V cir . Setting r c = R*^ (appropriate for a spherical 
geometry and a covering factor of unity), and using R* e ^ ~ 100 kpc, their model results give 
V C i r — 250 — 400 km s" 1 over 2 > z > 1, which is consistent with the observed kinematic 
ranges over 1.2 > z > 0.4 from Churchill & Vogt (2001). Using these values we obtain 
tinfaii ~ 0.25—0.4 Gyr. Thus, it would appear that a large degree of regulation is necessary, in 
this picture, to explain the lack of evolution in n(z) cr{z). Ionization from star formation and 
the disruption from instabilities and evaporation of smaller clouds should at least partially 
regulate the cooling. The result of these and/or other regulatory mechanisms must ensure 
that the minimum time span over which condensation occurs be equal to at least several 
times t infaU . 

Alternatively, lifetimes could indeed be much smaller than the time interval here studied, 
but a range of formation epochs conspire to maintain a roughly constant total cross section. 
While this scenario seems less likely than early formation epochs and long lifetimes through 
regulation, it is of note that although individual low-mass halos in ACDM simulations do 
themselves evolve, the low-mass end of the total halo mass function shows little evolution in 
redshift (Reed et al. 2003). 

The evolutionary situation for disks may be more complex. Simulations suggest that, 
to overcome "angular momentum catastrophe" , local disks were not formed until z < 1 (Mo, 
Mao, & White 1998; Weil, Eke, & Efstathiou 1998) and have since undergone much evolution 
(Mao, Mo, & White, 1998). Disks were likely plentiful at higher z, but smaller in size for a 
given V C i r , and largely dissipated or destroyed in mergers. Driver et al. (1998) find a deficit 



-20- 



of z > 2 spirals in the HDF suggesting that this marks the onset of their formation. Disks 
certainly contribute at least partially to the absorber population, as suggested by Charlton 
& Churchill (1998) and confirmed for specific cases by Steidel et al. (2001). However, even 
if their contribution to 0.3 A < Wq 2796 < 2.0 A Mg II absorption systems is significant, the 
lack of evolution remains difficult to explain. 



4.3.3. Evolution of Ultra-Strong Systems (W£ 279e >2A) 

We detect evolution in Wq 2796 Mg II systems, in the sense that the strongest lines 
evolve away, with the evolution being stronger for increasingly stronger lines and at redshifts 
z < 1. For example, we find that the total cross-section for absorption of Mg II systems 
with W A2796 > 2.0 A (« 10% of all systems with W^ 2796 > 0.3 A) decreases by 45l^% from 
z = 1.8 to 0.6. The total cross-section for absorption of systems with Wq 2796 > 3.5 A (ss 1% 
of all systems with W$ 2796 > 0.3 A) decreases by 60^% from z = 1.9 to 0.8. Note that for 
the Wq 2796 values considered in this work, Wq 2796 is correlated with the number of kinematic 
subsystems (clouds) comprising the absorber and thus velocity dispersion. 

A speciously simple picture for the evolution of Mg II A2796 lines considers that velocity 
dispersions scale with halo mass, and halo masses grow with decreasing redshift. In this 
picture, with the approximation that the cross-section for a given dispersion scales with mass 
and the mass spectrum evolves according to the formalism of Sheth and Tormen (Sheth & 
Tormen 1999), the cross-sections should increase with decreasing redshift by approximately 
an order of magnitude per unit redshift for 2 > z > 0. This is clearly ruled out by the data. 
This can be understood on physical grounds if the low-ionization gas giving rise to Mg II 
absorption is in clouds bound in galaxy halos, since the individual clouds are not expected 
to be virialized on group scales. A single Mg II system sampling a group dispersion would 
require the chance alignment of several virialized member galaxy halos along the line of sight. 

However, C2000 show kinematic profiles of four z > 1 and one z ~ system with 
Wq 2796 > 2 A, each of which exhibit two strong "dominant subsystem" components, and 
suggest the possibility that the strongest systems have a connection to galaxy pairs. For 
example, they find a 25% chance that a random sightline through the Galaxy would also 
intercept the LMC, and a 5% chance that it would intercept the SMC. This possible con- 
nection can be investigated by imaging Wq 2796 > 2 A systems: our sample contains 77 such 
lines with z < 1 and 20 with z < 0.6. 

Absorption dN/dz is the product of the number density of absorbers times the average 
individual cross-section for absorption. Most galactic halos are already in place by z = 1, 
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so except for mergers of galaxies separated by v > 500 km s^ 1 , the evolution in dN/dz is 
primarily an evolution of the gas absorption cross-section in individual halos. Evolution in 
metallicity or metagalactic ionization cannot explain the cross-section evolution, since they 
would tend to increase the number of individual enriched, low ionization clouds along a line 
of sight through a galaxy. The W / A2796 -dependent decrease in cross section must therefore be 
indicative of an evolution in the kinematic properties of a subset of the galaxy population, 
from intermediate to low redshift. If interactions, galaxy pairs and superwinds/superbubbles 
represent a significant fraction of the ultra-strong Mg II absorbers, then the decrease in the 
number of interactions/pairs from high redshift (e.g., in the HDF) and the decrease in the 
space density of superwinds (due to the decrease in the global star formation rate at z < 1) 
may in part account for the evolution in dN/ dz. Although these issues and those discussed 
in §4.3.2 likely hold clues to the nature of this evolution, a precise description will require 
knowledge of the physical nature of the clouds that comprise the different ranges of Wq 2796 
in Mg II absorption systems. 



5. Conclusions 

We have identified over 1,300 Mg II absorption systems with W^ 2796 > 0.3 A over the 
redshift range 0.367 < z < 2.269 in the SDSS EDR QSO spectra. The size of the sample 
is such that statistical errors are comparable to systematic effects and biases. We used 
simulations to improve our line-finding algorithm until the systematics could no longer be 
isolated from the noise at a level where we could further improve our parameterizations. 
Using the combined redshift and Wq 2796 data, we offered a new redshift parameterization 
for the distribution of systems. In conclusion, we have shown: 

1. The rest equivalent width distribution for intervening Mg II A2796 absorption lines 
detected in the SDSS EDR QSO spectra with W^ 2796 > 0.3 A is very well described by 
an exponential, with N* = 1.187 ±0.052 and W* = 0.702 ± 0.017 A. Power law param- 
eterizations drastically over-predict the number of strong lines, and our exponential 
under-predicts previously reported values for the number density of weaker lines. 

2. When compared with the number density of Wq 2796 < 0.3 A lines from other studies, 
our results show that neither an exponential nor a power law accurately represents the 
full range of Wq 2796 . Simulations of our catalog support this finding. We propose a 
combination of two exponential distributions, where the weak lines are described by 
the parameters N^ k 1.7 and W* k ~ 0.1 A and the moderate and strong lines by the 
parameters N* tr 0.9 and W* tr w 0.8 A. 
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3. The rest equivalent width distribution steepens with decreasing redshift, with W* 
decreasing from 0.80 ± 0.04 A at z = 1.6 to 0.59 ± 0.02 A at z = 0.7. 

4. For lines with Wq 2796 < 2 A, there is no significant evolution detected in dN/dz. 

5. For lines with W^ 2796 > 2 A, evolution is detected in dN/dz, with a decrease from the 
no-evolution prediction of ~ 45% from z = 1.8 to 0.6. The evolution is stronger for 
stronger lines and redshifts z < 1. 

6. The number density of Mg II absorption lines with Wq 2796 > 0.3 A is well parameterized 
by = e "TO, with Ar*(z) = 1.001 ± 0.132 (1 + z )^±o.no and = 

0.443 ± 0.032 (1 + ^)0-634±0.097 ^ 

Several lines of evidence suggest that the clouds giving rise to Mg II absorption comprise 
multiple physically-distinct populations. The results from high resolution work and the 
apparent transition in dN/dW^ 2796 are perhaps the strongest evidences, but the apparent 
contribution of enriched Lya forest lines, the inability of simple physical models to reproduce 
in full the kinematic data and the menagerie of galaxy types, luminosities, environments and 
impact parameters that contribute to the absorber galaxy population are also consistent 
with this picture. 

The situation will be made even clearer with the analysis of the full SDSS database. 
With over an order of magnitude more data available, there will be enough high signal-to- 
noise ratio spectra for studies to reach lines weaker than 0.3 A providing coverage across the 
ON/OWq 2796 transition in a single survey. Although the large number of systems will require 
large, improved Monte Carlo simulations, the improved statistics will allow for finer analysis 
of the strength-dependent redshift evolution. Additionally, work currently in progress will 
extend our analysis down to redshifts z ~ 0.2 and W^ 2796 > 0.1 (Nestor et al., in preparation). 
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Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics 
(MPA), New Mexico State University, University of Pittsburgh, Princeton University, the 
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A. Parameterization 

In order to parameterize d 2 N/dz dWo, we write 

d 2 N N*(z) Af* . a _ , 
= — e w z = ( 1 + z) p e v^y l+z > ( AD 

such that N*(z) = Af* (1 + z) a and W*(z) = W* (1 + zf '. The likelihood of the data set is 
then 

n(\ + /2 .)"-/3 e -^(i+^)- /3 
i 1 ^^; g 2 

i f f(l + z)«-Pe-^ 1+z ^ g(W X2796 ,z)dzdW ' 

The parameters that maximize the likelihood are a = 0.226 ± 0.170, f3 = 0.634 ± 0.097, and 
W* = 0.443 ± 0.032 A, which were determined with the aid of the MINUIT 4 minimization 
software. The maximum likelihood values are highly correlated, and the errors include the 
effects of the correlations. The resulting normalization is Af* = 1.001 ± 0.132. 

The full covariance matrix is 

0.0291 -0.00990 0.00322 

-0.00990 0.00931 -0.00292 j . (A3) 

0.00322 -0.00292 0.00104 

which should be used when calculating the uncertainty in d 2 N/dz dW - Af* was determined 
from the fit so that J J gfg^ 9(Wq 2796 , z) dzdW equaled the total number of lines in the 
survey. The full covariance array was used to calculate a^*. No contribution from the 
uncertainty in the number of lines was used, as it was estimated with a jackknife method 
to be small (< 3%). Since the uncertainty in Af* is derived from the uncertainty in the 
other parameters, it should not be considered when using equation Al. For example, for 







2 

a aW* 
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2 
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a W*a 


Cyy/3 





4 © CERN, Geneva 1994-1998 
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dN/dz = Af* (1 + z) a e~w T ( 1 + 2 ) the uncertainty is given by 

fd^Y 2 fdf^Y 2 f d^Y 
°dN/dz=<'L \ -^- I + < J ?/3 ( I +^v*w* I ^ 1 + ( A4 ) 

[w) [afc) +2a ^ t#j (^J ' 

where = f ln(l + *), ^ = f ^ (1 + ln(l + z), and |f = f (1 + *)"/». 
All uncertainties are statistical errors only. Possible systematics are discussed in §3.4. 
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Table 1. dN/dz Comparison 



w \2im 




This Work 




SS92 


Range (A) 


(z a bs) 


dN/dz 


(Zabs) 


dN/dz 


> 0.3 


1.11 


0.783 ± 0.033 


1.12 


0.97 ± 0.10 


> 0.6 


1.12 


0.489 ± 0.015 


1.17 


0.52 ± 0.07 


>1.0 


1.14 


0.278 ± 0.010 


1.31 


0.27 ± 0.05 
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Table 2. Summary of Data Parameterizations' 



N* 



W* (A) 



Full Sample: 



Data 
1.187 ±0.052 



0.702 ±0.017 



0.366 < z < 0.871 
0.871 < z < 1.311 
1.311 < z < 2.269 



1.216 ±0.124 
1.171 ±0.083 
1.267 ±0.092 



0.585 ±0.024 
0.741 ±0.032 
0.804 ±0.034 



Monte Carlo Simulations 



"weak" -phase: 
"strong" -phase: 



1.71 ±0.02 
0.932 ±0.011 



0.072 ±0.001 
0.771 ±0.014 



Redshift Dependence 



1.001 ± 0.132 (1 + z) 



0.226±0.170 



0.443 ± 0.032 (1 ± z) 



0.634±0.097 



a Maximum likelihood fits of the form: ON/OWq 2796 = w*. 
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The distribution of rest equivalent widths, Wq 2796 , for Mg II systems found in the 



Fig. 1. 

survey with W^ 2796 > 0.3 A. 
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Fig. 3. — Histogram of the Wq distribution. The solid line is a maximum likelihood fit of 
the form dN/dW^ 2796 = ^e~w* with W* = 0.702 ± 0.017 A and N* = 1.187 ± 0.052. The 
open triangles are from SS92. The dashed line is their best fit power law. The x symbols 
represent data from CRCV99. The dotted line is their power law fit to their binned data plus 
the SS92 data, excluding the highest-^ 2796 SS92 bin. The inset shows the W^ 2796 < 1.0 A 
region in more detail. The power laws greatly over-predict the incidence of strong (Wq 2796 > 
2 A) systems, while our exponential under-predicts the incidence of Wq 2796 < 0.3 A systems. 
This suggests a transition in the dN/dW^ 2796 distribution around W^ 2796 ~ 0.3 A, possibly 
indicative of multiple distinct populations. 
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Fig. 4. — Monte Carlo results. The histograms represent the data. The points are from 
CRCV99. The left panel shows the input exponential distributions. The shallower expo- 
nential (W* = 0.771 A and N* = 0.932) was chosen such that the simulated output distri- 
bution had W* ut = W* data for W A2796 > 0.5 A. The steeper exponential (W* = 0.072 A and 
iV* = 1.71) was then chosen to minimize the \ 2 computed by comparing the total simulated 
binned data to the real binned data. The solid line is the sum of the two exponential fits. 
The right panel shows the simulated output (solid line) over-plotted on the data. The dotted 
lines represent the square root of the variance from different random number seeds. 
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Fig. 5. — The redshift evolution of dN/dW^ 2796 . The solid lines are maximum likelihood 



fits of the form dN /dW^ = w 1 . Left: 0.366 < z < 0.871, with maximum likelihood 
value W* = 0.585 ± 0.024 A and corresponding N* = 1.216 ± 0.124. Center: 0.871 < z < 
1.311, with maximum likelihood value W* = 0.741±0.032 A and corresponding N* = 1.171± 
0.083. Right: 1.311 < z < 2.269, with maximum likelihood value W* = 0.804 ± 0.034 A and 
corresponding N* = 1.267 ± 0.092. Note the steepening of the slope with lower redshifts. 
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Fig. 6. — The redshift evolution of W*. The horizontal bars represent the bin sizes. The 
open circles are from the data, while the squares represent the Monte Carlo input values (see 
text, §3.1). The points are offset slightly in z abs for clarity. The solid line is the prediction 
for W* for a fit of the form W*(z) oc (1 + z) 13 described in §3.3. 
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Fig. 7. — The distribution of absorption redshifts for Mg II systems with Wq 2796 > 0.3 A 
found in the survey. 
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Redshift 



Fig. 8. — The total number of sightlines with sufficient signal to noise ratio to detect lines 
with W^ 2796 > W^ in as a function of redshift, for W^ in = 2.0 A, 1.0 A, 0.6 A, and 0.3 A. 
The conspicuous features at z > 1.5 are due to poor night sky line subtractions in many of 
the spectra. The depression near z— 1.1 is due to the dichroic. 
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Fig. 9. — Number density evolution of Mg II absorbers for Wq 2796 > W mm samples. 
The dashed lines are the no-evolution curves for a WMAP cosmology, (Qm,^a, ho) = 
(0.3,0.7,0.7). The curves are normalized to minimize the x 2 to the binned data. The 
W™« = 0.3 A, 0.6 A, 1.0 A, and 1.5 A samples have \ 2 values that are consistent with no 
evolution. The W^ in = 2.0 A, 2.5 A, and 3.0 A samples are inconsistent with the no evolution 
curves at > 3a, while the PV™ n = 3.5 A sample is inconsistent at ~ 2a. The dotted boxes 
represent the results of the Monte Carlo simulation. The widths correspond to the bin sizes 
and the heights to the ±lo~ values. The W m ' m = 3.0 A and 3.5 A samples were not large 
enough to permit meaningful Monte Carlo simulations. 
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Fig. 10. — Number density evolution of Mg II absorbers for ranges of Wq 2796 . The dashed 
lines are the no-evolution curves for a WMAP cosmology, (Om, &a, h ) = (0.3, 0.7, 0.7). The 
curves are normalized to minimize the x 2 to the binned data. The 0.6 A - 1.0 A, 1.0 A - 1.5 A, 
1.5 A - 2.0 A, 2.0 A - 2.5 A and 2.5 A - 3.0 A samples have x 2 values that are consistent with 
no evolution. The no-evolution curve for the 0.3 A - 0.6 A sample is ruled out at ~ 2.5a and 
for the 3.0 A - 3.5 A and > 3.5 A samples at ~ 2.0a. The dotted boxes represent the results 
of the Monte Carlo simulation. The widths correspond to the bin sizes and the heights to 
the ±la values. The 3.0 A - 3.5 A and greater than 3.5 A samples were not large enough to 
permit meaningful Monte Carlo simulations. 
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Fig. 11. — The cumulative distribution of dN/dz (solid lines) for the Wq 2796 > W™ in 
samples, compared with the cumulative no-evolution curves (dashed lines) for a WMAP 
cosmology (f2 M , f2 A , /i ) = (0.3,0.7,0.7). The dotted curves represent the la levels. The 
no evolution curves under-predict dN/dz at low redshift for the W™ m = 0.3 A sample, and 
over-predict dN/dz at low redshift for the larger W™ m samples. The data are inconsistent 
with the NECs at more than 3a for all but the W™ m = 0.6 A sample. 
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Fig. 12. — The cumulative distribution of dN/dz (solid lines) for the W^ 2796 ranges, com- 
pared with the cumulative no-evolution curves (dashed lines) for a WMAP cosmology 
(£l M ,il A ,h ) = (0.3,0.7,0.7). The dotted curves represent the la levels. The no evolu- 
tion curve under-predicts dN/dz at low redshift for the W^ 2796 = 0.3 A - 0.6 A and W^ 2796 = 
0.6 A - 1.0 A samples, and over-predicts dN/dz at low redshift for the larger Wq 2796 ranges. 
However, the data rule out the NECs at > 3a only for the Wq 2796 > 2.0 A samples. 
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Fig. 13. — dN/dz curves calculated from our parameterization (§3.3). The dN/dz data are 
consistent with the fit for all ranges shown. Left: The cumulative samples. Right: The 
non-cumulative samples. 
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Redshift Redshift 

Fig. 14. — The cumulative distribution of dN/dz plotted against the cumulative dN/dz 
curves calculated from our parameterization (§3.3). The dN/dz data are consistent with the 
fit for all ranges shown. 
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Fig. 15. — Left: the distribution of W 7 "^ 2796 /]^ 2803 doublet ratio for Mg II systems with 
Wq 2796 > 0.3 A. The dashed lines mark the limits of 1.0 for completely saturated systems 
and 2.0 for completely unsaturated systems. Values above and below these limits are due to 
noise. Our sample is dominated by saturated systems. Right: Wq 2803 versus Wq 2796 . 
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Fig. 16. — Wq 2796 /Wq 2803 doublet ratio as a function of redshift for Mg II systems with 
W^ 2796 > 0.3 A. There is no detected redshift evolution in the doublet ratio. 
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Fig. 17. — Cumulative distribution of approximate velocity dispersions for <J ve i jP > 90 km 
s -1 . The lower and middle redshift bins differ in K-S probability at the ~ la level. 
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Fig. 18.— Mg II - Fe II ratio. Left: The distribution of Mg II A2796 - Fe II A2600 line 
ratios. The largest bin also contains all systems with ratios greater than 8.0, which includes 
non-detections of Fe II A2600. The mean is (W 7 ^ 2796 /!^ 2600 ) = 1.42 and ratios above « 4 are 
dominated by values with significance less than 3a. Right: Wq 2600 versus Wq 2796 . Negative 
yj/A2600 va i ues are non-detections with scatter below zero. 
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Fig. 19.— Mg II - Mg I ratio. Left: The distribution of Mg II A2796 - Mg I A2852 line 
ratios. The largest bin also contains all systems with ratios greater than 25, which includes 
non-detections of Mg I A2852. The mean is (W^ A2796 /W^ A2852 ) = 4.14 and ratios above w 8 are 
dominated by values with significance less than 3a. Right: Wq 2852 versus W A2796 . Negative 
yj/A2852 va j ues are non-detections with scatter below zero. 
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Fig. 20. — Predicted ON/OWq 2796 for Lya forest/single-cloud Mg II . The solid line is the 
maximum likelihood fit to our data with W^ 2796 > 0.5 A. The dotted line is the approximation 
for single-cloud Mg II dN/ dW^ 2796 using data from CRCV99 and a normalization from Rigby, 
Charlton, & Churchill (2002). The dot-dash curve is the sum of the two distributions. The 
points are from CRCV99. 
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Fig. 21. — Cross-section estimates cr(L*) = nR* e ^. The solid lines represent the L min = 
0.05 L* result. The dashed-lines are for a small value of L min = 0.001 L* and the dotted- 
lines a large value L min = 0.25 L*. We use the redshift parameterization of the K-band LF 
from the MUNICS data set, a scaling relation R e ff/R* e ^ = (L /L*) , and dN/dz from our 
parameterization (§3.3). 



